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Abstract. TeV y-rays from the BL Lac object 1ES 1959+650 have been measured during the years 2000 and 2001 with a 
significance of 5.2 cr at a value of 5.3 % of the Crab flux and in May 2002 during strong outbursts with > 23 <x at a flux level 
of up to 2.2 Crab, making 1ES 1959+650 the TeV Blazar with the third best event statistics. The deep observation of 197.4 h 
has been performed with the HEGRA stereoscopic system of 5 imaging atmospheric Cherenkov telescopes (IACT system). 
1ES 1959+650 is located at a redshift of z =0.047, providing an intermediate distance between the nearby Blazars Mkn 421 and 
Mkn 501, and the much more distant object H1426+428. This makes 1ES 1959+650 an important member of the class of TeV 
Blazars in view of the absorption of TeV photons by the diffuse extragalactic background radiation (DEBRA). The differential 
energy spectrum of 1ES 1959+650 during the flares can be fitted by a power law with a spectral index of 2.83 ± 0. 14 stat ± 0.08 sys 
or by a power law with an exponential cut-off at (4.2*9 - ? ± 0.9 sys ) TeV and a spectral index of 1.83 ± 0.15 stat ± 0.08 sys . The 
low state differential energy spectrum obtained with lower statistics can be described by a pure power law with a spectral index 
of 3.18±0.17 stat ±0.08 sys . 

Key words, y-rays: observations - BL Lacertae objects: individual: 1ES 1959+650 



1. Introduction 

Active Galactic Nuclei (AGN) are known to be sources of ex- 
tragalactic TeV y-radiation. Except for a recent report of a sig- 
nificant TeV excess from the radio galaxy M 87 observed with 
the HEGRA Cherenkov telescopes (Aharonian et al. 2003b I all 
TeV AGN detected so far are of the BL Lac class. In these 
objects the very high energy photons are believed to originate 
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(possibly due to inverse Compton scattering) in the relativis- 
tic jets oriented at small angles to the observer's line of sight. 
To these TeV y-ray emitters belong the well studied BL Lac 
objects Mkn 421 (z = 0.030) and Mkn 501 (z = 0.034). 

Recently, the much more distant BL Lac H1426+428 
(z = 0.129) has been firmly established as an emitter 
of TeV y-radiation by the VERITAS dHoran et al. 20021 
[ Petty et al. 2002), HEGRA JAharonian et al. 20021 
Uharonian et al. 2003a> . and CAT ( |Djannati-Atai et al. 2002) 
collaborations. Possible further TeV detections of BL Lac 
objects have been reported by different collaborations for 
lES2344+514(z = 0.044) dTataiiese et al. 19981, the BL Lac 
PKS2 155-304 (z = 0.117) dChadwick et al. 19991 and the 



2 



Aharonian et al.: Detection of TeV Gamma-Rays from the BL Lac 1ES 1959+650 



very distant object 3C 66A (z = 0.444) ( |Neshporet al. 19 98 1, 
but no confirmation by any other group was stated so far. 

The object 1ES 1959+650 Elvis etal. IffHl with a red- 
shift of z = 0.047 was classified as a BL Lac in 1993 
using a specially developed radio/optical/X-ray technique 
( Schachte r et al. 1993i . A first detection at TeV energies was 
reported by the Utah Seven Telescope Array collaboration for 
the 1998 observational season ( Nishi yama et al. 1999| l. An ex- 
cess with a statistical significance of 3.9 cr above 600GeV 
was obtained after 57 h. There was no independent confirma- 
tion until the year 2001, when the HEGRA collaboration re- 
ported a detection of 1ES 1959+650 with a significance of 
> 4cr (|G otting et al. 2001) , resp. > 5 cr with a larger data set 
IHornsetal. 20021 1. 

In May 2002, 1ES 1959+650 underwent a strong TeV 
outburst, detected by the VERITAS ( Holder et al. 2002> . 
HEGRA JHorns et al. 20021 IKranich et al~2002l and CAT 
dVorobiov et al. 2002b collaborations. Strong variations with 
flux levels of up to 3 times the Crab flux have been measured. 

In this Letter the results of extensive observations 
of 1ES 1959+65 using the HEGRA IACT system 
JDaum et al. 19971 during the years 2000/2001 and 2002 
are reported. Astrophysical conclusions concerning the nature 
of this fourth established TeV Blazar regarding possible 
influence by the DEBRA absorption are briefly discussed. 

2. Observations and results of analysis 

For the following analysis, a total of 18. 4h of data from the 
year 2000, 89.4h from 2001, and 89.6h from 2002 were used 
taken with the HEGRA IACT system. The dates of the HEGRA 
observation periods are listed in Tabled The mean zenith an- 
gle of all observations was 38.8°, resulting in a mean energy 
threshold (defined as the peak detection rate for y-showers) of 
1 .4 TeV for a Crab-like spectrum ( Konopelko et al. 1999 1. The 
1ES 1959+650 data were mainly taken with the complete 5- 
telescope setup of the IACT system. A few data runs were taken 
with only 4 telescopes due to technical reasons. Individual runs 
were accepted, if well-defined quality criteria were fulfilled 
(see, e. g. Aharoni an et al. 2003b> . About 8 % of all data were 
rejected due to this selection. 

All observations of 1ES 1959+650 were carried out in the 
so-called wobble mode allowing for simultaneous estimation of 
the background ("OFF") rate induced by charged cosmic rays 
( Aha ronian et al. 19971 1. This analysis uses for the signal search 
the so-called ring segment background model as explained in 
lAharonian et al. 20""0 3b providing a small ratio of ON to OFF- 
source solid angle areas a = A ON /A OFF (resulting in improved 
background statistics). For the spectral analysis a set of 7 inde- 
pendent background regions is used ( Aharo nian et al. 2001> . 

The general shower reconstruction and the event selection 
cuts for the image analysis have been described in previous 
publications (e. g. lAharonian et al. 199% . For the signal search, 
the stereo air shower direction reconstruction algorithm #3 
( Hofman n et al. 199 9 1 and a tight shape cut (parameter mscw 
< 1.1) (Konopelko et al. 1999 1 for an effective y-hadron sepa- 
ration have been applied. Because of the relatively large zenith 
angles more usable shower images per event are recorded corn- 



Table 1. Individual periods of 1ES 1959+650 observations 
with the HEGRA IACT system. Listed are the total observa- 
tion times and mean zenith angles {ff). Typically, each night 
comprises about 1 h of observation time. 



Date 


Year 


Obs. Time 


\ u / 






Thl 

L U J 


L J 


July 28 - August 2 


2000 


8.2 


40.0 


August 28 — September 6 


2000 


10.2 


40.7 


May 1 8 — June 1 


2001 


13.0 


37.4 


June 22 - June 28 


2001 


15.5 


38.1 


July 1 1 - July 22 


2001 


8.5 


38.1 


August 9 - August 25 


2001 


17.5 


37.2 


September 6 - September 19 


2001 


15.0 


37.3 


October 1 1 - October 21 


2001 


19.9 


40.1 


May 18 -May 21 


2002 


10.0 


41.8 


June 3 - June 18 


2002 


34.4 


38.5 


July 1 - July 18 


2002 


18.7 


39.5 


July 31 - August 16 


2002 


20.2 


38.4 


September 3 - September 1 1 


2002 


6.3 


38.6 


Sum 


197.4 


38.8 



Table 2. HEGRA IACT system direction reconstruction algo- 
rithms and cuts used for the 1ES 1959+650 signal search and 
for the spectral analysis, respectively. The resulting numbers 
for 4 selected subsets are given: "2000/2001" including all data 
used for the first 5.2 cr detection, "May 2002" only using data 
from the nights May 18/19 and May 19/20, 2002 (see Figure^. 
The "low state" resp. "high state" data sets are used for the 
spectral analysis (see text and Figure[3}. 





signal 


spectrum 


stereo algorithm 


#3 


#3 


number of images per event 


>3 


>3 


shape cut on mscw 


< 1.1 


< 1.1 


angular distance cut O 2 


< 0.012 deg 2 


< 0.014 deg 2 


a = AonMoff 


0.0758 


0.1429 




2000/2001 


"low state" 


Tobs 


93.8h 


149.6 h 


A'ON 


354 


585 


A'oFF 


3441 


1960 


A^ y-candidates 


93 ± 19 


305 ± 25 


significance (cr) 


5.2 


14.6 




May 2002 


"high state" 


Tobs 


3.6h 


8.5 h 


A'ON 


179 


275 


A'oFF 


131 


141 


A^ y-candidates 


169 ± 13 


255 ± 17 


significance (cr) 


23.4 


25.5 



pared to observations close to the zenith. Therefore, a minimum 
number of 3 images per event is required for the whole recon- 
struction chain rejecting all 2-telescope events, which only pro- 
vide a poor resolution. The optimum angular cut was derived 
using y-ray events from the Crab nebula on the basis of a nearly 
contemporaneous data set at similar zenith angles. The proce- 
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dure of spectral evaluation - with an energy resolution of < 
12 % for a single event - is discussed in Aharo nian et al. 2003al 
The energy spectrum is derived from a background subtracted 
photon count spectrum using an effective area on an event by 
event basis, depending upon the number of active telescopes, 
the reconstructed energy, the zenith angle and the assumed 
spectral shape. In order to take into account the slightly vary- 
ing detector performance, the effective areas are determined on 
a monthly basis along with the mirror reflectivities, photocath- 
ode efficiencies and conversion factors from digitized photo- 
multiplier signals to photoelectrons ( Puhlho fer et al. 2003I I. 

The relevant direction reconstruction algorithms, event se- 
lection cuts, resulting event numbers, and significances for the 
individual HEGRA 1ES 1959+650 data sets used for the sig- 
nal search and the spectral analysis are summarized in Table|2 
Figure ^ shows the event distribution for the ON-source and 
the OFF-source regions from the May 2002 flare and the year 
2000/2001 quiescent state data sets, respectively, as a func- 
tion of the squared angular distance to the source position. 
The statistical significance of the excess from the direction of 
1ES 1959+650 amounts to 5.2 cr in the year 2000/2001 obser- 
vations, calculated using formula 17 in Li & Ma (19 83) The 
evidence for TeV y-rays from 1ES 1959+650 during this qui- 
escent state at a value of (5.3 + l.l stat )% of the flux of the 
Crab nebula (Aharonia n et al. 2000 ) already confirms the ten- 
tative detection of 1ES 1959+650 by the Utah Seven Telescope 
Array. The HEGRA detection during the strong outburst in 
May 2002 is obvious at a very high significance of greater than 
23 cr. The diurnal integral flux levels above 2 TeV observed in 
the year 2002 are shown in Figure |2] (Table 4 containing the 
complete lightcurve for the years 2000 - 2002 is available at 
the CDS, see note to the title). The strong outburst in May 2002 
is followed by a post-flare low state in June. During July, the 
source has exhibited more activity again followed by a period 
of diminishing TeV flux in August and September 2002. 

The spectral investigation of the 1ES 1959+650 "high 
state" was performed using the data of the 6 nights in May 
and July 2002 with the object's integral flux being larger than 
1 Crab above 2 TeV (see Table [3j. Within the limited statis- 
tics, the single nights of these 8.5 h of observations show the 
same spectral shape. A fit of a pure power law dN/dE = 
No x (E /lTeVy a to the differential energy spectrum results 
in N = (7.4 ± 1.3 stat + 0.9 sys ) x 10" 11 phot, cirT 2 s^T/eV^ 1 
and a = 2.83 + 0.14 stat ± 0.08 sys with ^ (d.o.f.) = 1.9(6). 
A fit of a power law with an exponential cut-off dN/dE = 
No x (E/lTeV)-°> x exp(-E/E c ) with N = (5.6 ± 0.9 stat + 
O^xlO-Hphot.cm-V^eV- 1 ,^ = 1.83±0.15 stat +0.08 sys , 
E c = ( 4 -2 tLa, ± 0.9 sys ) TeV, ^ (d.o.f.) = 1.7(5), is also an 
adequate description of these data (see Fig. [3}. 

Figure [3] also shows the spectral energy distribution of 
the 1ES 1959+650 "low state". Nearly 150h of observations 
during nights with 1ES 1959+650 emitting TeV y-rays at a 
flux level below 0.5 Crab have been used for this purpose. 
The spectrum can be described by a pure power law with 
N = (7.8 + 1.5 stat + 1.0 sys ) x lO-^phot.cirrVTeV- 1 and 
a = 3.18 + 0.17 8ta , + 0.08 sys with X l A (d.o.f.) = 0.22(3). A 
power law with an exponential cut-off and a spectral index fixed 
at a = 1.8 fits the data with N = (6.0 + 1.4 stat + 0.8 sys ) x 
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Fig. 1. Number of events vs. squared angular distance O 2 to the 
position of 1ES 1959+650, observed on May 18/19 and 19/20, 
2002, with the HEGRA IACT system. The dots show the ON- 
source events, the histogram gives the background estimate de- 
termined from a ring segment as explained in the text. The sta- 
tistical error for the background estimate is much smaller than 
the error of the ON-source distribution. Upper right: Similar 
distributions for the year 2000/2001 observations yielding a 
significance of 5.2 cr. 

Table 3. HEGRA IACT system event statistics and values of 
the differential photon flux of 1ES 1959+650 in the high state 
for the individual energy bins as shown in Figure [3] 



Energy 


Won 


NoFF 


Sign. 


diff. Flux 




(TeV) 






(<T) 


(phot, cirr 2 s" 1 TeV" 


') 


1.52 


42 


21 


10.0 


(2.5 ± 0.5) x 10" 


ii 


2.02 


44 


43 


8.6 


(8.5 ± 1.6) x 10" 


12 


2.68 


44 


25 


10.0 


(4.0 ± 0.7) x 10" 


12 


3.55 


54 


12 


12.9 


(3.0 ± 0.4) x 10" 


12 


4.70 


34 


9 


10.0 


(1.3 ± 0.3) x 10" 


12 


6.23 


16 


3 


7.1 


(4.5 ± 1.2) x 10" 


13 


8.25 


6 


4 


3.5 


(1.1 ±0.5) x 10" 


13 


10.94 


5 


1 


4.0 


(0.7 ± 0.3) x 10" 


13 



10~ 12 phot.crrT 2 



s 

2 (d.o.f.) = 0.65 (3) 



TeV" 1 and E c = (2.7 



+o.6 +06 ") TeV 



3. Summary and Conclusions 

The BL Lac object 1ES 1959+650 has been observed with the 
HEGRA IACT system for a total of 197.4 h. The object has 
been detected at TeV energies during its low flux state in the 
years 2000 and 2001 at a value of (5.3 + l.l stat )% of the Crab 
flux and during flaring states in 2002, reaching a flux level as 
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HEGRA-1ES1 959+650 



(1 Crab) 



May 2002 i June 2002 



52500 52520 
modified Julian date (MJD) 

i August 2002 i Sept. 2002 



10h 
23.0 a 



34h 
6.4 



19 h 
23.6 a 



19 h 
16.2 a 



3h 
2.5 a 



Fig. 2. Diurnal integral flux values resp. upper limits 
(99 % c. 1.) observed above 2 TeV vs. the modified Julian date 
(MJD) for the year 2002 HEGRA IACT system data. Strong 
flux variations are visible in May and July. The level of the flux 
of the Crab nebula is indicated by the horizontal dashed line. 
The time gaps between the individual observation periods are 
due to the full moon periods when no observations were possi- 
ble. 



high as 2.2 Crab. Extensive HEGRA observations in 2002 show 
strong changes in the absolute flux level. 

The time averaged energy spectrum of the high state obser- 
vations is well described by a pure power law with spectral in- 
dex a - 2.83 +0.14 slal + 0.08 sys or by a power law with an expo- 
nential cut-off at E c = (4.2+° ® ^ ± .9 sys ) TeV. Within the syn- 
chrotron self-Compton (SSC) model the TeV y-ray emission of 
a BL Lac type object corresponds to the Inverse Compton (IC) 
component. |Costamante & Ghisellini (2002)| have predicted 
the IC flux above 1 TeV in a flaring state of 1ES 1959+650 
to be 1.74 x 10" 11 phot, cm 2 s 1 using a phenomenologi- 
cal parametrization of the spectral energy distribution by 
|Fossati et al. (1998)1 Taking into account the strong variations 
of 1ES 1959+650 in 2002 the observed flux values can be eas- 
ily accommodated. 

Due to the absorption of TeV photons by the pair produc- 
tion process y TcV + Tdebra — * <? + + e~ with the DEBRA the ob- 
served spectrum differs from the source spectrum. In the range 
from 1 to several TeV however, the spectral shape remains 
nearly unchanged since the optical depth in this energy range 
is rather constant for most of the DEBRA candidate spectra in 
the relevant wavelength range of 1 to 10 fim (Aharonian 2001 1. 
Using a "model of choice" prediction of the DEBRA spectral 
energy density (lAharonian 200 It . the attenuation coefficients 
for the HEGRA 1ES 1959+650 high state data points have 
been calculated and applied to unfold the measured spectrum as 
shown in Figure[3] The cut-off energy is only marginally shifted 



to (4.5 ^ g stat + 0.9 sys ) TeV. The low event statistics at the high- 
est energies observed around 10 TeV and above (Tab.[3j do not 
allow one to see the strong deformation of the spectral shape, 
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Energy [TeV] 

Fig. 3. Time averaged spectral energy distributions (SED) of 
1ES 1959+650 as measured with the HEGRA IACT system. 
The filled circles show the spectrum during 6 high state nights 
in 2002, while the filled squares show the SED for the combina- 
tion of all low state nights (< 0.5 Crab). The results of fits of a 
power law with an exponential cut-off (solid lines) resp. a pure 
power law (dashed lines) are indicated. The open circles show 
the high state spectrum corrected for the DEBRA absorption 
(see text). 

which is expected from the energy dependence of the attenua- 
tion coefficients (e. g. see Figure 9 in Aharo nian 200 It . 

The fact that 1ES 1959+650 is about 1.5 times more dis- 
tant than the nearby TeV Blazars Mkn 421 and Mkn 501 and 
bridges the distance gap to H1426+428 makes it a very impor- 
tant object for future DEBRA probing in the optical to near- 
infrared region. 
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